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A generalized multibaker map with periodic boundary conditions is shown to model boundary- 
driven transport, when the driving is applied by a "perturbation" of the dynamics localised in a 
macroscopically small region. In this case there are sustained density gradients in the steady state. 
A non-uniform stationary temperature profile can be maintained by incorporating a heat source 
• into the dynamics, which deviates from the one of a bulk system only in a (macroscopically small) 

' localized region such that a heat (or entropy) flux can enter an attached thermostat only in that 

, region. For these settings the relation between the average phase-space contraction, the entropy flux 

04 ' to the thermostat and irreversible entropy production is clarified for stationary and non-stationary 

(— I ' states. In addition, thermoelectric cross-effects are described by a multibaker chain consisting of 

, two parts with different transport properties, modelling a junction between two metals. 

I. INTRODUCTION 

' There is a recent interest in modelling transport processes by simple dynamical systems with chaotic dynamics. 



X 



One class of models, actually inspired by non-equilibrium molecular dynamics (NEMD) simulations, describes sys- 
tems driven by external fields with a spatially uniform dynamics subjected to periodic boundary conditions [yj-[7|- 
^ ' Another approach concentrates on systems driven from the boundaries, which lead to steady states with sustained 
. gradients of the thermodynamic fields |^,^ . For a comparatively simple, but as far as their transport properties are 
concerned, generic class of dynamical systems, the multibakers, pO|-pO| we show that both mechanisms of driving can 
simultaneously be worked out. This leads to an improved understanding of the relation between the approaches. In 
7-H ■ the former approach transport is driven by a field acting uniformly in the full system, while in the latter case the 
^ , driving is concentrated to a microscopic region in space. From this point of view, boundary-driven transport is closely 
^^O ■ analogous to transport in a dynamical system with periodic boundary conditions, which is driven out of equilibrium 

' by a "perturbation" of the dynamics, localised in a macroscopically small region. 
^"^ In all models for transport, as emphasized by Nicolis and coworkers |^,^, a quantity of central interest is the heat 

flux, or equivalently the entropy flux, from the system into its environment. A central aim of modelling transport by 
, dynamical systems is to identify settings, which are consistent with the thermodynamic entropy balance 

' d5 dpS diS . . 

.9: dF"^^^' 

Ch ' i.e., with the statement that the temporal change of the thermodynamic entropy S can be decomposed into two 
\ contributions, called the external and internal change of the entropy, respectively. This integral form can be rewritten 
into a local balance equation when the two terms on the right hand side correspond to integrals of local densities. In 
that case, the time derivative of the entropy density s appears as 

-'^ i 9ts = $ + cr^"---) (2) 

where $ and cr^"''') represent the densities of the entropy flux and the rate of irreversible entropy production, respec- 
tively. In the bulk of typical macroscopic systems the entropy flux can be written as the divergence of the entropy 
current j 

<i>--V/^\ (3) 

reflecting the fact that no heat can be taken out from the system locally |2|]. On the other hand, this form has to 
be generalized at positions where there is a heat current flowing into an attached thermostat, and in cases where the 
entropy current is not differentiable, like for istance across interfaces between different materials. In those cases the 
entropy flux is not a full divergence, and it need not even be defined as a density. Rather the flux should then be 
written as 

-Vj(") -f (4) 
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where j^^^ is still the entropy current flowing in the system, but accounts for the heat taken out in the form of 

a direct flow into the surroundings, which acts then as a thermostat. 

In the present paper, we shall put special emphasis on the role of the entropy flux and on exploring under 

which conditions it can vanish. We find conditions on how to model the entropy balance for thermodynamic bulk 
systems, and for macroscopic systems, which are subjected to thermostatting by either a localized sink for the entropy 
or a spatially uniform coupling to a thermostat. 

The role of the thermodynamic entropy of dynamical systems is played by the coarse-grained Gibbs entropy, whose 
usefulness in understanding irreversibility from the point of view of dynamical systems is by now thoroughly explained 
in the literatu re l22|- 27 , 1^-|20| . ( For stochastically perturbed dynamical systems where noise generates a kind of coarse 
graining, see P4y25[ ). The bulk dynamics is represented by a multibaker model driving two fields, the density g and 
the the kinetic energy per particle T, with a local source density q for the latter A connection with macroscopic 

transport equations is aimed at in a suitable defined continuum limit (the macroscopic limit), where the field T will 
be interpreted as a temperature, based on the experience that this quantity is closely related to the average kinetic 
energy per particle. 

We shall consider a sequence of periodic models of increasing complexity. Model I corresponds to a homogeneous 
isothermal system described by a thermostatting algorithm. In this model no entropy current is defined — its entire 
entropy flux stems from a By allowing a spatial resolution of the isothermal system (Model II), a non- vanishing 

— Vj^''^ term appears in the transient behaviour, but remains unchanged. It is the only contribution to the flux 
in a steady state. Model III is still isothermal but with a locally deviating dynamics in one of the multibaker cells 
representing a boundary. The bulk dynamics can then be chosen so that (^) holds in the bulk, and all the heat taken 
out is concentrated in the boundary with a ^ there. In model IV we allow for temperature changes and local 
heat sources. By taking q locally deviating from that of the bulk in one cell, we find a steady temperature profile with 
a break at the boundary. The q distribution can then be chosen such that again (^ holds in the bulk. The heat source 
in the boundary is however singular. It corresponds exactly to the one which follows from thermodynamics. Finally, 
we consider a multibaker chain joined together from two subchains with different material properties. This models a 
junction between two metals so that one can observe thermodynamic cross effects, like the Peltier and Seebeck effects, 
very much in the same arrangements as in classical experiments. 

This paper is organized as follows. In Section II the local dynamics of the considered multibaker model is defined, 
and its local entropy balance is worked out. In Section III, Model I - Model IV are treated, which share periodic 
boundary conditions and represent thermodynamic settings of increasing complexity. Section IV is devoted to cross- 
effects. We conclude with a short discussion in Section V. 
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FIG. 1. Graphical illustration of the action of the multibaker map on cell m. The letters L, S and R are inserted to 
visualize this action. Iteration of the rule after every time unit r defines the time evolution. More details about the action of 
the mapping and the symbols needed for its definitions are given in the text. The symbols gi and Ti indicated on the margins 
show the average values of the fields on the cells and on its neighbors i — m i 1. 

II. LOCAL TRANSPORT AND THERMODYNAMIC RELATIONS FOR MULTIBAKERS 



In this section we describe the local dynamics of a cell of a multibaker map modelling a system with particle 
pO| , p^ 14 -1^ and heat transport [^Q. We work out its density and kinetic-energy dynamics, and present general 
relations for the entropy changes. The effect of boundary conditions will be considered in subsequent sections for a 
few models with progressively richer thermodynamics. 

The phase space {x,p) of the multibaker map consists of cells labelled by the index m (Fig. The division of 
the X axis into cells corresponds to a partitioning of the configuration space into regions, sufficiently large to allow 
to characterize the state inside the cell by thermodynamic variables and small enough to neglect variation of these 
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variables on the length scale of the cells (local equilibrium approximation). Every cell has a width a and height b = 1. 
The coordinates of individual particles in the cell are given as a position variable x, and a momentum- like variable 
p. We are interested in the dynamics of two dimensionless fields, the phase-space density g{x,p), and a field T{x,p) 
characterizing the local kinetic energy per particle. After each time unit r, every cell is divided into three columns 
(Fig. ^) with respective widths al„i, aSm and arm- (Note that Im + Sm + Tm = ^ for any to.) The right (left) column 
of width ar„i {aim) is uniformly squeezed and stretched into a strip of width a and of height Im+i (''m-i), which is 
mapped to the right (left) neighbouring cell. The middle one preserves its area and remains in cell m. Note that 
the map is one-to-one on its domain. It globally preserves the phase-space volume, but it can nevertheless locally 
expand or contract the phase-space volume. In Ref. it was argued that only the choice of the contraction 

factors given here can be consistent with thermodynamics (in fact, one can find an analogous formulation with a fully 
area preserving dynamics at the expense of a spatial variation of the volume of the cells of the multibaker; cf. p9|). 

The field T is advected by the particle dynamics, and — in order to mimic a local heating of the system — it is also 
multiplied by a factor (1-1- Tq) depending on the averages characterizing the local currents and the thermodynamic 
state. By this a mean- field-like coupling of the motion of the particles in and around of a given cell is introduced. In 
general, the width of the columns may depend on the variables characterizing the thermodynamic state in the vicinity 
of the cell, so that they vary in time and space. This is indicated by the explicit dependence of the parameters on 
the cell index. Iteration of these rules defines the time evolution of the system. The (x,p) dynamics generates ever 
refining structures in the distributions g{x,p) and T(x,p). For simplicity, we take the fields initially constant in each 
cell: g{x,p) = T{x,p) = . 

A. Dynamics of the particle density and the particle current 

After one step of iteration, the fields will be piecewise constant on the strips defined in Fig. |l|. Due to the conservation 
of particles, the phase-space density takes the respective values 

/ ^m — 1 / / /|-\ 

Qm.r ~ ~i Pm-1, Qm.s ~ Qrn, Qm.l ~ Qm+1- \p) 

(The prime will always indicate quantities evaluated after one time step.) The contraction factors rm-i/lm and 
Im+i/rm contribute to the (weighted) local phase-space contraction Um of cell to: 



Pm-lfm-l In — h Qm+llm+l In 



(6) 



After one time step, the average density g'^ in cell to is determined by its initial density gm and by the initial densities 
gm±i of the neighbouring cells. Multiplying the strip densities (^) with the widths of the respective strips, adding 
them up and dividing the sum by the width a of the cell, one obtains the average (or the coarse-grained) density after 
the iteration 

dm — Sm.Qm + ?'m-lgm-l + Im+lQm+l- (J) 

The coarse-grained density evolves according to this master equation, which can be rearranged to obtain the discrete 
conservation law of the density 



6m 9m jm jm — 1 



(8) 



Here 



jm = -{rmQm - lm+l6m + l) (9) 



T 



is the discrete particle current flowing through the right boundary of cell to. 

B. The kinetic-energy dynamics and the energy current 

According to the T dynamics described above, the updated values ^, T^^, ; for T on the respective strips 
R, S, L of cell to contain a source term characterized by a local strength qm- 
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TL^s = Tm [1 + rg™] , (10) 
T!m,i = r,„+i [1 + rq,„] . 

This strength is yet undetermined. It depends on the physical setting of thermostatting to be modelled and on the 
average of g and T values in the cells and in its neighbours. 

The (a;,p) dynamics also drives the field T, i.e., after one iteration the kinetic-energy density of cell m takes the 
value 

T„+i](l + Tg„). (11) 
This equation can be rearranged as a discrete balance equation for the time evolution of qT: 

n T' n T a i^^^^ - l^^^^ 

= QmTm-r~ ' 12) 

T 1 + a 

where jm"^'* = Tmjm — {a^ Itu+i/t) Qm+i (Tm+i — Tm)/a is a corresponding discrete energy current. Note that the 
r.h.s of ( |l2|) is not a full divergence, in accordance with the fact that the kinetic energy is not a conserved quantity. 
In an isothermal system where there is no kinetic energy dynamics, no source can be present (qm = 0). 



C. Gibbs entropy and the coarse-grained entropy 

In this study we are interested in both the temporal evolution of the exact fields g{x,p) and T{x,p), and in the 
evolution of their respective cell averages gm and T^. The former densities characterize the microscopic time evolution, 
while the averages describe the local thermodynamic state in spatially small regions. Both levels of description admit 
entropy functionals, which are commonly denoted as Gibbs and coarse-grained entropy. 

The Gibbs entropy S**^*^^ is related to the detailed knowledge of the system. It is taken with respect to the exact 
densities g{x,p) and T{x,p). In a given cell it is defined as 



= - [ dxdp g{x,p) In f T {x , pY 

J over cell m \ 



(13) 



Here g*T'' plays the role of a local T-dependent reference density with a constant reference density g* and 7 an as 
yet undetermined constant. 

The coarse-grained entropy S„i has a similar form, but it is based on the averaged values in the considered cell: 

S„, = -ag,nln(^^T-,-'y (14) 

As mentioned above, throughout the paper we only consider initial distributions, which are uniform in every cell 
(cf. |l^,|l^,Q for more general choices). As a consequence, initially Sm = Sm\ and after one time step the entropies 
become 



and 



Srngni In T„^J^ + r'm-lfi'm-l In ^ T„i,r^ + ^m+l(?m+l In ^ T,n,l^ 



(15) 



S'^^^ag'„Jn[^T:-^]. (16) 
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D. Entropy balance 



The coarse-grained entropy fulfills a local entropy balance in direct analogy to the one in irreversible thermody- 

(G) 

namics. To derive this equation one identifies at any given time the difference Sm — Sm as the information on the 
microscopic state of the system, which cannot be resolved in the coarse-grained description. The temporal change of 
this lack of information is then identified with the irreversible entropy production AiSm, and the change {Sm — Sm ) 
of the Gibbs entropy with the entropy fiux Ai,Sm- Thus, 



S^n Sn 



Af^Sri 



AiSn 



(17) 



which is a discrete analog of (|^). 

The form of the entropy production is [cf. (|l5|) and (p^)]: 



\S' - S- 



(G)'i 



[5„ 



S 



(G)i 



dm In 



^m m 



m^s 



Qm-lTm-l In 



~m,r m.r 



QmTj 



/-7 



+ Qm+llm + 1 In 



9m.lT, 



I 



QmT, 



m^s 



(18) 



(G) 

where we used that Sm — Sm vanishes due to the particular choice of initial conditions. 
The entropy flux becomes 



In ( ^T„^ 



ftnln 



Qm—l'^n 



In 



Qm Tm.s 



+ Qm+llm+1 In 



which can be split into a divergence of an entropy current and a flux into the thermostat 

AeSm 



Jm Jm~l 



dm.l '^m.l 



(19) 



(20) 



with 



$(th) = _i 

m — 

T 



rp -7 

e'm In + Tm-lQm^l In 



'^mQm '^m — lQm—1 



a 



Qm'^rt 



+ Un+lQm+l In 



'; 7^-771-7 



(21a) 



(21b) 



Note that ( |2C| ) is a discrete counterpart of and jm ^ and <^^m^ are the discrete entropy current and entropy flux 
to the thermostat, respectively. 



E. The macroscopic limit 



The projection of the multibaker dynamics on the x axis corresponds to a biased random walk with some diffusion 
coefficient and drift. The drift has to be present if we want to model nonequilibrium systems subjected to electric 
fields and/or temperature gradients. The requirement of consistency with an advection diffusion equation in the large 
system and long time limit, when the cell size is much smaller than the system size, and the time unit is much shorter 
than the macroscopic relaxation time, leads |12,nJ,na,n^ to the scaling relation: 



fm 5" I 

I -^f 

'^m — o \ 



1 



1 - 



2D 



(22a) 
(22b) 
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Here we allow for a location dependence of the drift Vm but assume the diffusion coefficient D to be spatially constant. 
The continuum limit of the multibaker dynamics, which is taken with these constraints, is called the macroscopic 
limit. Formally, it corresponds to taking a,r — > while keeping D fixed, and am, Qm, Tm, Vm and qm finite so that 
they approach a macroscopic position coordinate x, and smooth functions g{x), T(x), v{x) and q{x), respectively. 
After taking this limit we call g{x) the density and T(x) temperature distribution in the system. The macroscopic 
limit of all the local relations given in (0)~(|l^) can be worked out explicitly [[I2yi4|jl5| , [l7|] . In the following this limit 
will be indicated by an arrow '— >' . Here we only mention that the system's equation of state turns out ]l7| , p^ to be 
that of classical ideal gas with as its heat capacity (measured in units of Boltzmann's constant). 



q=0 T=const 



111 I r I 

FIG. 2. Isothermal single-cell periodic baker map. The right boundary is identified with the left one. 

III. PERIODIC MODELS 
A. Model I: Isothermal single-cell multibaker 

We start by discussing the simplest conceivable model for describing a macroscopic transport process. A particle 
current induced by an external field in an isothermal environment described by a single baker cell subjected to periodic 
boundary conditions. The right boundary of the cell is identified with its left boundary and the mapping is from the 
cell onto itself. Because of driving r ^ I, and thermostatting is applied via the appearance of the contraction rates 
l/r and r/l in order to reach a steady state. This mapping propagates the coordinates of a large number of particles, 
which do not interact, i.e., they all are mapped by the same mapping. Clearly, this system does not admit a spatial 
resolution of the densities characterizing the transport process. Its local and global behaviour coincides, so that the 
subscript m of the densities can be descarded in this case. 

For the single-cell multibaker the master equation (0) predicts g' = g = g. This implies that the model is describes 
tranport in a steady state with the average density g. The particle curent j = {a/T)(r — l)g = v g is constant in space 
and time, and the entropy production (|l§| ) becomes (T = const.). 



which has the macroscopic limit: 



a^-'-) = -g-. (24) 



In a similar way, the entropy flux (gO|, pib| ) has the macroscopic form: 

M ^ $(th) ^ _ V 

ar D 

As expected in a steady state, AiS and AeS add up to zero. More interestingly, however, these contributions to the 
change of entropy are also directly proportional to the local phase-space contraction (^, which reduces to gv'^/D in 
the macroscopic limit. 
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FIG. 3. Isothermal multibaker chain of TV + 1 cells. The right boundary of cell TV is identified with the left boundary of cell 



B. Model II: Isothermal multibaker chain for field driven transport 

In the single-cell multibaker the coarse-grained density cannot evolve in time, since no spatial variations are resolved. 
For this reason, one also cannot distuinguish between the local entropy balance as described by (^), (^) and the balance 
(|l|) for the full macroscopic system. In order to address these points, we generalize the previous setting by considering 
a multibaker chain of iV -I- 1 cells, with spatially constant driving (i.e., Vm = f and Im = and periodic boundary 
conditions: cell iV -t- 1 and cell are identified. From (0) we obtain for the particle current 



3n 



- [{r - l)Qn 
T 



'm+1 



which has the macroscopic limit: 



(26) 



(27) 



This current varies along the chain as long as the cell densitites evolve in time. The asymptotic state, however, is 
formed by a spatially uniform density distribution with the average density q. 
For constant temperature T the irreversible entropy production (|l^) becomes 



I 1 



'•Qm-l In 



iQm+l In 



±1 

, r 



(28) 



Using (0) and (p^), a lengthy but straightforward calculation for a, r — > shows that this local form is consistent 
with thermodynamics since in a general non-steady state it approaches 



{gv - DdxO) 
qD 



2L 

qD 



(29) 



in the macroscopic limit. In a similar way [ p^ , p^ , the entropy flux ( p^ can be calculated, which in the macroscopic 
limit has the form (M), with the entropy current 



j(^) =-j [l + ln(^,/e*)] 



(30) 



and the entropy flux 



{th) ^J^± 
D 



(31) 



transfered directly to the environment. This local expression expresses that every cell is coupled to the thermostat. 

Due to the additional spatial resolution (as compared to Model I) the local and global features of the entropy 
balance can be different. It is worth considering the global entropy production, defined as [cf. (28)] 



N 

E 

m=0 



A,S„ 



N 



— 



l)\n 



N 
rn—Q 



[-gi'„ln g'^ + g„i\n g„ 



(32) 
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Here, the logarithm of the ratios of densities of neighbouring cells in (^8[) drop out in the sum due to the periodic 
boundary conditions. The global entropy production then takes the macroscopic form 



N 
m=0 



djS 

dt 



■ TV 

E 



2 ^ 



m—0 



D 



dS_ 

dt ' 



where M — ag{N + 1) is the total number of particles in the multibaker chain, and S = 
The global form of the entropy flux is [cf. (|l9|)] 



(33) 



Sm is the total entropy. 



N 



ApS'r, 



m=0 



N 



(r - In ( 



N 



m=0 



which takes the macroscopic limit 



m 



deS 

dt 



-4- 



(34) 



(35) 



This shows that the local and the global entropy balances are markedly different. Locally, the entropy current depends 
only on the local current and the local density. There is only an indirect influence of the drift velocity v through its 
contribution to the current and its influence on the density profile. In contrast, the macroscopic flux only depends on 
the total number of particles in the system and on the drift velocity. It is constant in time since it neither depends 
on the current nor on the density profile, which in general both evolve in time. Interestingly, also in this more 
general setting the negative of the global entropy flux equals the (total) phase space contraction at any time. This 
is in full harmony with the result obtained for the entropy flux in noisy dynamical systems by Nicolis and Daems 
1^,^. In contrast, the total irreversible entropy production is in general no longer directly related to the phase- 
space contraction. Its local form exactly amounts to Joule's heating / gD, and globally it picks up an additional 
contribution characterizing the time-evolution of the macroscopic states. Therefore, the often cited relation between 
the global entropy production and the phase-space contraction only holds in steady states, where the modulus of the 
entropy flux coincides with the rate of entropy production. 



q=0 T=const 



-//- 
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FIG. 4. Isothermal multibaker chain with driving in cell 0. There is no drift in the bulk. 



C. Model III: Isothermal multibaker chain with driving localized to a single cell 

A generalization of the previous case is a system modelling boundary-driven transport. Here, we consider a slightly 
different setting with periodic boundary conditions. Away from the point of driving the macroscopic properties 
correspond to that of boundary driven transport, but due to the periodic boundary conditions the model has a 
simpler structure as a dynamical system, which makes the analogy with our previous models more transparent. 

We consider a multibaker chain where cell has a behaviour different from the rest of the chain in the sense that 
its drift velocity is different from that of the bulk where r = I, (i.e., v = 0) but the diffusion coefficient is the same. 
This implies Tq — r ^ I — Iq ^ (cf. (^2|)). The driving in cell zero has the tendeny to generate an accumulation of 
particles right to it and a slowly decreasing density distribution on its left (we assume Vq > 0, rg > r). In a steady 
state, this leads to a linear density profile 

gm = g + ( — ■;; rn + 1 I 6g for m = 1...N, (36) 
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whose increment 5q is uniquely determined by tq (or vq). In cell zero we find go = g. A substitution into the master 
equation for m = 1 or leads to 

^.= .-^(^-1). (37) 

Again the average density g is related to the number of particles J\f in the system via ag{N + 1) = JV. By taking into 
account that according to ( |22| ) rg = r[l + avo/{2D)], we obtain that 6g is indeed proportional to vq: 

SgN+l SgN 

Vq = D ; = £> —. (38) 

a g a ag^ 



The steady state current is: 



-Vo T^Qi^S ^Sg vog 
2 a a N + 1 



The diffusion current of cells N and is much stronger than otherwise, due to the sudden jump in the densities, but 
in the steady state the surplus is exactly compensated by the local drift currents. Thus the particle transport in the 
steady states of Model II and Model III are equivalent provided the drift v in the former coincides with DSg/{ag) in 
the latter. As far as the steady state transport is concerned, it does not matter, whether one applies a small uniform 
field leading to a spatially uniform drift w or a large one vq — {N + l)v in a single cell only. 

All thermodynamic relations of relevance can be worked out not only for the steady state (^6|), but for general 
non-steady states. The local forms of the entropy production and the entropy flux in the bulk are special cases of 
(p9|)-(|3l|). Since r = Z in the bulk, there is no entropy flux flowing into the thermostat — for m = 2, ...N — 1. 

On the other hand, the entropy flux in cell and its neighbours does contain a part which cannot be written as a 
divergence. Due to (^lb|), the average density of the flux flowing into the thermostat turns out to be: 



^ Vq gi- gN Vq gl + gjV Vq .^qx 

^ 2 a AD 2 AD' ^ ' 

In the last approximation we have assumed that avo ^ 21?. 

In the steady state gjv ~ gi = ^(A^ ^ l)'^g = ^(A^ ~ l)aj/-D, go — ''^^^^"'^ — g, and thus 

D ' ^ ' 

In the macroscopic limit the quantity Sg/a approaches the gradient —dxg in the bulk, g and A/" stay constant and 
thus Vq in ( ^ ) is proportional to 1/a. The driving is singularly strong and so is the entropy flux density into the 
thermostat. By integrating, however, over the volume of cell zero and its neighbours we obtain the total entropy 
flux into the thermostat —av^j/D = —j'^M/g^D. It coincides with the macroscopic limit of the total entropy flux 
X]m=o ^eSm/T siucc the integral of dxj^"'^ vanishes in a periodic system, i.e.. 



dt gD g 



(42) 



This result is equivalent to the steady state version of ( |35| ) expressed by the current. 

In this model there is no need for taking out heat along the bulk, thermostatting is active in cell and its neighbours 
only. It extracts exactly the same entropy flux there as the full entropy flux of Model II in the steady state g] . Thus 
the models refer to two different realizations of thermostatting the transport process. Model II should be viewed 
as e.g. a wire which is kept at constant temperature by removing the heat due to dissipation everywhere along its 
length — Model III is closer related to a thermally isolated system, where heat is tranported to the "boundary" , from 
where the system is driven. For the multibaker this takes place in the special cell m — Q. Boundary driven transport 
typically leads, however, to non-uniform temperature profiles. A full treatment of such transport processes should be 
based therefore on a multibaker chain with kinetic energy dynamics. 
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FIG. 5. Non-isothermal multibaker chain with thermostatting localized to cell 0. The local heat source go of that cell is 
diflerent from that in the bulk, which is q* = —vj/D. 



D. Model IV: Multibaker chain with thermostatting locahzed to a single cell 

In the realm of classical thermodynamics heat is transported to the boundaries of the system. If, however, there is 
a break in the temperature profile, a jump in the heat current occurs and heat is taken out at this point. 

In order to model such a situation, we consider a multibaker chain with kinetic-energy dynamics. The general 

relation (|l^) ^and a calculation similar to the one leading to (|2^) yields in the macroscopic limit for the entropy 

production Q 

+ (43) 

where j = gv — DdxQ is the particle current, and A = ^qD is the heat conductivity of the model. From (|l9|)-(21b) we 
obtain 

^^^''^^TOT-^ (44) 
as the entropy flux let directly into the thermostat, and 

= -A^ + f (45) 
as the entropy current with 11 the bulk Peltier coefficient 

f^-(>..»^). 

Note that there is always a possibility to 'close' the system locally in the sense that the source term q — q* ~ vj/X is 
chosen such that vanishes. 

We consider a periodic chain with fixed transition probabilities (r„i = r, Im = rn = 0.1, -.jN). The local heating 
sources are assumed to be constant in the bulk: q^ = q* for m — 1, N which differs from the source qo of cell 0. In 
the steady state we find then a constant particle density along the chain. Inside the bulk, the kinetic-energy equation 
( |rT| ) implies for the steady temperature distribution 

T„, - [(1 - r - l)T„-, + rTra^i + ff„+i](l + rq*). (47) 

With periodic boundary conitions (T,„=o = ^Oi a-nd Tm=N+i — Tq) this equation has the following general solution: 

m /- W + 1 -V 

T„ = -[^(1^ (0 " [mHN + 1) - + (1) ' sin(6.n)| (48) 

where 

cos6 = Vwfl-- -V (49) 
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The solution (g8|) has a break in cell zero, in the sense that the left and right derivatives are different. Only at the 
end of the system is the entropy flux not a full divergence. Applying the kinetic-energy equation ([ll| ) to cell zero in 
a steady state, where the density is constant, we find that: 



^ 1 r(ro-rAr) + ^(ro-ri) 

^° T (1 - r - l)To + tTn + ITi ■ 



In the macroscopic limit 



(-0) 



(+0) 



(50) 



(51) 



This implies that the source density is singular but the total source Qo — ago inside cell 0, is finite. 

It is worth comparing this with the thermodynamic treatment of the same problem. If there is a jump in the 
entropy current, in order to have finite entropy flux density $ in each point, it is unavoidable to allow for a $ which 
is not a full divergence: 



-a,/^) + -I. 



(52) 



The form of $ one obtains by integrating (52) around the point where the jump in the derivative appears {x — 0): 



The smoothness of <I> implies that for e — > 



/: 



^dx. 



(+0) 

(-0) ■ 



(53) 



(54) 



The multibaker result ( |44| ) implies that if q is singular as in cell 0, then = 'ygqo- We then immediately see that 



<I>Q is the analog of $. Indeed, by inserting the expression (45) for j^'^^ we have: 



T\ 



(-O)J 



(55) 



which, on account of A = ^qD, exactly corresponds to (pO|). 

We have shown, that the thermodynamic evaluation and the macroscopic limit of go lead to the same result. 
Physically this means, that by a proper choice of the source terms even the singularity in the entropy flux can be 
described in full harmony with thermodynamics, and the flux let to flow in the thermostat is exactly the amount of 
heat what is taken out also in the thermodynamic description if a break appears. 



-II- 



-II- 



material 

-//- 



material + 



in I r I 



■II- 



■II- 



11 1 I r II I I r t 1^1 I r I 1 I I r I 



■II- 



■II- 



I r I I r* I 



■II- 



m= -M' ••• -2-10 1 ■■■ M 

FIG. 6. Two long multibaker chains, representing materials — and +, joined together at the junction between cells —1 and 
0. The leads are in ceUs -M' and M. 



IV. CROSS EFFECTS IN A MULTIBAKER MODEL 



Thermodynamic cross effects, which probe the (off-diagonal) Onsager coefficients, are difficult to observe in ho- 
mogenous systems. When two materials are put into contact, however, they play a dominant role in understanding 
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the heat and entropy currents. In order to mimick such phenomena, we consider two multibaker chains containing 
the cells m = —M' • • • — 1 and m = • • • M, respectively, M, M' >> 1 which are brought into contact at m = 
(cf. Fig. H). Now, the parameters s~ , r~ and s^, r+ in the two parts are different, and for generality we 
will also assume that the constant reference densities g*'^ are different. These differences will represent the different 
thermodynamic and transport properties of the materials. The difference in r and I gives rise to different drifts (con- 
ductivities) and diffusion coefficients, and the one in the reference density might be thought of reflecting for instance 
a different effective mass of the electrons. 

As in the previous subsection, the dynamics of this multibaker chain drives a density and a kinetic-energy field. 
In order to simplify the structure of the steady state density profiles, we restrict to the case r+/Z+ = /l~ . This 
choice is motivated by a physical interpretation of r/l. After all, the macroscopic limit of r/l is 1 -I- av/ D, and v/D 
is proportional to the external (electric) field, such that the requirement expresses that the external field should be 
the same in both materials. In the remainder of this section we discuss the transport in this model in two different 
settings: (i) a constant (non-vanishing) particle current and constant temperature; (ii) vanishing particle current and 
an isolated system which is only thermostatted at the "junction" m = and at the two "leads" m = —M' and 
m = M, respectively. Setting (i) allows us to discuss the Peltier effect and setting (ii) is used for the Seebeck effect. 

Before turning to these specific settings, however, we discuss the steady state profile of the (particle) density in 
general One does not expect noticable gradients in the electron density in either material, so that we fix them to the 
constant values g~ and g'^ , leading to the spatially uniform current (cf. (^) 

J = -{r+-l+)9+ = -{r- -ng-. 

T T 

In order to have the same current also across the junction, one has to require l~ g'' ~ g'^ in adition. Together with 
the fact that v/D is fixed for the whole system, this implies that there is a constant amount of Joule's heating vj/D 
per unit length of the system, which either has to be transferred to a local thermostat (cf. (^)), or leads to a local 
heating, i.e., enforces non-uniform temperature profiles. 



A. The Peltier effect 



The requirement of a constant temperature in the setting of the Peltier effect requires the use of a thermostatted 
dynamics, q = 0. In that case the Joule heating is transferred to the thermostat. Away from the junction, this leads 
to the flux = vj/D. In the entropy balance the difference in the materials shows up only in the entropy currents. 
In view of (Esh , they become different in the two parts of the multibaker 



J 



-J 



1 + In 



-J 



implying that at the junction (i.e., between cells m — —1 and m = 0) an additional heat flux, the Peltier heat, is 
directed to the thermostat. It is characterized by the difference of the entropy currents 



J 



J 



-j In 



g+/g 
g- /g*' 



J 



In- 



F 



In- 



en(+/-) 



T 



where 11'+/ ' defines he mutual Peltier coefficient of the two materials. It characterizes the amount of Peltier heat 
produced per unit electric current, and is the difference of the material Peltier coefficients [cf. (B^)] 



n(+/-) = n+ - n- 



(56) 



as also found in thermodynamics. 



B. The Seebeck effect 



The Seebeck effect is observed in a thermally isolated system, where the junction is kept at a temperature Tj 
different from the temperature Ti prescribed at the leads, i.e., for the multibaker we demand T-m' = Tm = Ti and 
r_i = Tq = Tj. This setup corresponds to a non-uniform temperature field, and, due to this, also to gradients in 
the electro-chemical potential Because of the difference in the material properties, these gradients can add up to 
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a net potential drop between the leads, even if both leads are kept at the same temperature and there is no particle 
current. This follows immediately from the formal definition of the particle current in its discrete version: 



~ ~ ^2 



m+l 



(57) 



where Uei is the conductivity, e the electric charge, and a the Seebeck coefficient of the material. In the considered 
system, we then have for vanishing current 



ss -ea^ [T^M' - T^i) + M-i - Mo - ea+ {Tq - Tm) = e(Q!+ - a^) [Ti - Tj) + - ^o- (58) 

Here we have assumed the Seebeck coefficients to be approximately constant in the two materials. The macroscopic 
limit implies = /^Oj ^^nd we obtain for the mutual Seebeck coefficient of the two materials pSf 



(+/-) _ M-A/' - MAf + 

" e{Ti-T,) 



(59) 



It characterizes the strength of the potential drop H-m' ~ IJ-m between the leads induced by the temperature diffence 
T; — Tj between the leads and the junction. 

An expression for a'+Z"' can be determined for the multibaker by rewriting the expression (|2^) for the current in 
the form (^^. Taking immediately the macroscopic limit and observing that the electro-chemical potential can be 
split into a chemical part /^c and a part ecf) due to the external electric field E = —dx(f>, one obtains 



J 



Tel E 



[^g^J■c dxQ + dxt^c dxT + ea d^T] ^ vg- D d^g. 



(60) 



Here 



V 

D 

ea 



(7 el E 

eg 

f^el ^ 



(61a) 

(61b) 
(61c) 



By the first two equations we recover well-known relations from thermodynamics [ pi[ . Eq. (61c) provides us with a 
relation for the Seebeck coefficient. Since the equation of state of the "multibaker gas" is that of a classical ideal gas 



^it = (7 + i)r + rin 



one obtains 



v(+/-) - 



= e{a —a ) = e In 



g- 



g- 



.•1 



In- 



Ti 



where ( p6| ) was used in the last step. This comparison expresses the validity of the Onsager relation H^+/ ^ = a'+Z 
for this class of models. 



V. DISCUSSION 



In this paper we have described the local and global transport properties of multibakers with a density and an energy 
dynamics. This class of maps makes an analytical modelling of transport processes by a deterministic chaotic dynamics 
possible, and admits a macroscopic description consistent with various aspects of irreversible thermodynamics. The 
macroscopic description comprises the time evolution of the average density and the kinetic energy in small regions 
of the physical space (the cells of the multibaker). The former density is interpreted as the particle density, and the 
latter as a temperature field. The averages in the small regions are in the spirit of local thermodynamic equilibrium, 
and the continuum description of thermodynamics arises in a macroscopic limit where the spatial resolution of the 
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tranport process is small compared to the system size (or any other relevant macroscopic length), and where a discrete 
time-scale used in the definition of the dynamics is much smaller than macroscopic time scales. 

The relevant concept of entropy for multibakers is the Gibbs entropy defined with respect to the average density in 
the cells normalized by a temperature-dependent reference density. It is called the coarse-grained entropy. Based on an 
information-theoretic interpretation of the entropy, a local entropy balance can be derived, which in the macroscopic 
limit can be fully consistent with irreversible thermodynamics. This agreement holds provided that (i) a particular 
choice of local phase-space contraction and expansion rates is incorporated in the time evolution of the density, 
which we identified as a time-reversible evolution of the mapping in previous work (ii) the density in 

the entropy is normalized by a reference density with a power-law dependence on the average kinetic energy in the 
cell, and (iii) appropriate source terms are incorporated in the evolution equations of the kinetic-energy field. No 
meaningful macroscopic description can be found for multibakers with other choices of the phase-space contraction 
factors. Modification of the source terms leads to additional contributions in the local entropy balance, which are 
interpreted as local entropy fluxes into a thermostat. In particular, for vanishing source terms one can mimic a 
transport process in a system with a spatially uniform temperature, i.e., one obtains a setting reminiscent of NEMD 
simulations of transport processes. 

Once the connection between the deterministic dynamics of the multibaker and the corresponding local thermody- 
namic relations is established, one can apply it to discuss transport in different macroscopic settings. A number of 
models with periodic boundary conditions were discussed in order to shed light on the global entropy balance in such 
systems. We find that, up to a trivial factor, the average phase-space contraction amounts to the entropy flux to the 
environment. This supports an earlier heuristic argument of Ruelle |Q and others [^-||, who connected the phase- 
space contraction to the irreversible entropy production in a steady state. In contrast to the claims of some of the 
latter authors (cf. for instance p9|]) the connection between the irreversible entropy production and the phase-space 
contraction rate breaks down away from stationarity. In fact, the contraction rate is still connected to the entropy flux 
in that situation, but the flux is no longer related to the rate of irreversible entropy production. This was shown (a) 
for multibakers with a uniform thermostatting (Model II), i.e., for models reminiscent of NEMD algorithms, (b) for 
systems where the driving and thermostatting is applied in a macroscopically small region of the system (Model III), 
thus giving rise to sustained density gradients, and (c) systems with a uniform external field and localized thermostat- 
ting (Model IV). The former two models have constant temperature fields, while the latter one supports a temperature 
profile with a discontinuity in the first derivative at the positon of thermostatting. As expected from the existence 
of the local entropy balance, the results are fully consistent with the corresponding thermodynamic description of 
the transport process. They suggest an interesting conclusion on modelling transport in bulk systems by isothermal 
NEMD simulations: These methods are valid in an approximation where the considered volume is sufficiently small 
to neglect density and temperature gradients. In steady states, they are equivalent to models, where the currents are 
the same, but thermostatting is only applied at the boundaries of a macroscopic system. Since even state of the art 
simulations can hardly cope with more than 10^ particles, i.e., with integration volumes larger than about //m^, this 
approximation seems to be well-justified in numerical studies. On the other hand, this assumptions should be kept in 
mind when isothermal NEMD modelling is taken as basis of theoretical studies of transport processes (cf. for instance 

Mil)- 

To further demonstrate the use of multibakers with density and energy fields, we also discussed thermoelectric cross 
effects. The description of the transport properties requires in that case information on the equation of state, since the 
Seebeck effect is defined in terms of differences of chemical potentials. In previous work [[l7|jlq ] it was shown that the 
classical ideal-gas equation holds for multibakers. This is meaningful since the time evolution of the multibaker can 
be considered as the one of particles with phase-space coordinates {x,p), which only interact by a (weak) mean-field 
like coupling manifested in a dependence of the local parameters on the average densities. With this input the Peltier 
and Seebeck effect were modelled and the Onsager relation, connecting their respective transport coefficients, was 
derived. The validity of this relation for multibakers is not a trivial result. It heavily relies on the choices (i)-(iii) to 
find an entropy balance consistent with irreversible thermodynamics. 

Summarizing, we demonstrated that multibakers establish a straightforward modelling of various transport phe- 
nomena by deterministic, chaotic dynamics. They give insight in the general structure of such models by explicit 
analytical calculations. This was demonstrated by discussions of thermoelectric cross effects, and of the relation 
between the average phase-space contraction, entropy fluxes and the rate of irreversible entropy production. 
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